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In a recent work [H. Flayac, LA. Shelykh, D.D. Solnyshkov and G. Malpuech, Phys. Rev. B 81, 
045318 (2010)], we have analyzed the effect of the TE-TM splitting on the stability of the exciton- 
polariton vortex states. We considered classical vortex solutions having cylindrical symmetry and 
we found that the so-called half vortex states[Yu. G. Rubo, Phys. Rev. Lett. 99, 106401 (2007)] 
are not solutions of the stationary Gross-Pitaevskii equation. In their Comment [M. Toledo Solano, 
Yu.G. Rubo, Phys. Rev. B 82, 127301 (2010)], M. Toledo Solano and Yuri G. Rubo claim that 
this conclusion is misleading and pretend to demonstrate the existence of static half-vortices in an 
exciton-polariton condensate in the presence of TE-TM splitting. In this reply we explain why this 
assertion is not demonstrated satisfactorily. 
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First of all, we now agree with the authors of the 
Commenlii that our claim of the non-stability of half- 
vortices^ in the presence of TE-TM splitting is not sup- 
ported enough. Indeed, in our paper— we considered 
only solutions of the type ip± (r, <f>) — f± (r) exp (id± (</>)), 
where the amplitudes f± of the wave function depend 
only on r and phases 0± are linear with respect to <p. 
From this point of view the publication of a comment 
seems reasonable. 

On the other hand Solano and Rubo considered only 
asymptotic equations. They found nontrivial half-vortex 
asymptotic solutions and we agree with their result. In 
the other limit — > 0, the kinetic energy dominates all 
other terms, and a cylindrically symmetric half-vortex 
solution can be found. However, in the intermediate re- 
gion (close to the vortex core), all three contributions of 
the Hamiltonian (kinetic energy, interactions energy, and 
TE-TM splitting) are of the same order. The existence of 
the asymptotic solutions close to infinity and close to the 
ground state does not guarantee the existence of a half- 
vortex solution in the whole space. As we show below, 
the amplitude of these asymptotic solutions depends on 



both r, and (j). The general solution is therefore of the 
type ip± (r, (j)) = f± (r, tfi) exp (i6± (r, </>)). The analysis of 
the existence, or non-existence in the whole space of such 
type of half vortex is a very demanding task, which has 
not been performed neither by the authors of the Com- 
ment, nor by the authors of the reply. Qualitatively, we 
hardly believe that such vortices, with particles moving 
on non-circular trajectories could exist as their symmetry 
differs from that of the Hamiltonian and therefore some 
preferential directions appear. The authors of the Com- 
ment believe the opposite. At this stage, there is still no 
definitive answer given to the question of the stability of 
half vortices in the presence of TE-TM splitting. 

The asymptotic equation was obtained by the authors 
of the Comment by variation of the elastic energy of the 
vortex neglecting completely the radial dependence. We 
believe that this way important terms could be lost. The 
starting point should not be the elastic energy (which is 
not the total energy), but the Gross-Pitaevskii equations, 
which in the dimensionlcss form read (see Eqs. 11, 17 of 
Refi): 
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where according to the notations of the Comment we used </> f° r the polar angle and for simplicity neglected 
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the interactions of polaritons with opposite circular po- 
larizations, a2 = 0. We divided all the terms into three 
groups: those containing derivatives by <j> only (group 
1), those without derivatives (group 2) and those con- 
taining derivatives by r only (group 3). Considering the 
asymptotic case, one sees that the terms of the group 1 
are of the order of magnitude of r -2 , while the terms of 
the group 3 are of the order o(r~ 2 ) and thus can be ne- 
glected, as it was done by the authors of the Comment. 
As for the terms of the group 2, it can be easily shown 
that they are of the same order of magnitude as terms 
from the group 1. Indeed, using the approximative so- 
lution for the radial function of a vortex with a winding 
number 1, Eq. 23 of Refill one easily sees 
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Thus, the asymptotic equation reads: 



where: 
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If the solution is cylindrically symmetric, A± should be 
independent on <j). 

Now, using the ansatz of the authors of the Comment 
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one has 4 equations for the unknown functions 9, 77, A± 
(since both real and imaginary parts of the equation [3] 
should be equal to zero): 
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The main equations of Toledo Solano and Rubo (Eqs. 
6a and 6b of the Comment) can be obtained taking the 
sum and the difference of Eqs. |H] and [7] above, which 
ensures that the imaginary part of Eq[3]is zero. However, 
to make the real part being zero as well, one needs to 
satisfy additionally Eqs. 18191 as well. These equations, 
not written in the Comment, allow to deduce the phi 
dependence of \±. 

For solutions of the class found in Refi, A± arc- evi- 
dently constant and given by r\ = <f>, 9 = p<p, 9" = if = 0, 
?/ = 1,9' = p with p being an integer number and 
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On the other hand, if 9 and 77 are complicated non- 
linear functions of (j>, in order to satisfy Eqs l8l9l A± is 
cj> dependent, which means that the amplitudes of the 
vortex solutions at large but finite distance from the core 
region are non-cylindrical, cf> dependent: 
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where the dependence of X± on <j) is deduced from 
Eqs l8l9l As one clearly sees, in these solutions the ro- 
tational symmetry present in the initial Hamiltonian is 
broken and preferential directions appear. The peculiar 
choice of the author of the Comment is that the polariton 
pseudo-spin is aligned with the TE-TM effective field at 
<p = 0. But any other direction could be chosen equiv- 
alcntly. The existence of the solutions of this type in 
the whole space is not a priori excluded, but is neither 
guaranteed. 

Regarding the qualitative remark of Toledo Solano and 
Rubo that singularities in Gross-Pitaevskii equation can 
not appear or disappear we would like to mention that 
first of all, for a multi-component (e.g. spinor) conden- 
sate, in contrast with a scalar one, the conservation of 
current circulation and, consequently, vortex stability, is 
not in general guaranteed by topological considerations. 
In such a system a vortex can be removed by a con- 
tinuous transformations, which has been demonstrated 
experimentally^. Also it was already shown that vor- 
tices and vortex-antivortex pairs can appear in Gross- 
Pitaevskii equation from regular initial conditions, e.g. 
in a process of scattering of a BEC on localized defect^ 
or in the rotating condensates^. It was also shown that 
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TE-TM splitting can generate vortices of a winding num- 
ber 2 from regular initial distributions of the polariton 
condensate^. 

To conclude, we believe that the claimed stability of 
half-vortices under the effect of TE-TM splitting was 
not fully demonstrated. We disagree with the remark 
of the authors of the Comment that our article contains 
mathematical error as we found the exact vortex solu- 
tions possessing the symmetry of the Hamiltonian. The 
existence of half- vortices in the presence of TE-TM split- 
ting is not proven by the authors of the Comments, and 
this statement should not be claimed by them. On the 



other hand, some of the conclusions of Ref. are indeed 
not fully supported. Further works will be needed in or- 
der to give a definitive answer to this problem. Finally, 
one of the main conclusion of our work was that under 
the influence of the TE-TM field specific pairs of half- 
vortices become bounded to form the (— 1, +1) integer 
vortex solution (this was also remarked by the authors 
of the Comment). As a half- vortex will never live alone 
inside the microcavity one should reasonably expect that 
such a bounding phenomenon will occur and should leave 
the BKT transition temperature linked with integer vor- 
tices. 
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